Dissipative Landau— Zener Tunneling at Marginal Coupling 
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The Landau-Zener transition in a two level system can be suppressed or enhanced by coupling to an 
environment, depending on the temperature and the environment spectral function. We consider the 
marginal spectral function, when the dissipation effects are important for arbitrarily slow motion. 
Landau-Zener transition rate demonstrates a non-trivial dependence of the on the "bias", i. e., on 
the rate of the two energy levels relative motion. The Landau-Zener transition is fully suppressed 
for the values of the bias below a threshold bias set by the coupling strength. Above the threshold, 
the transition rate for zero temperature is found using the instanton method. At finite temperature, 
the Landau-Zener transition rate has a non-monotonic dependence on the coupling strength, being 
suppressed at the strong coupling. 
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Introduction. The Landau-Zener problem |^ deals with 
non-adiabatic transition in a two-state system under ex- 
ternal bias. Physical examples include current states 
in small metallic loops 1^, spin tunneling in magnetic 
molecules Andreev states in SNS junctions slow 
atomic and molecular collisions electron transfer in 
biomolecules Q. Real- world two-state systems are typi- 
cally coupled to the environment, and this coupling may 
lead to different physical effects. This coupling provides 
dissipation and slows the tunneling down. On the other 
hand, thermal noise in the environment may cause non- 
adiabatic transition. The competition between these two 
effects makes the problem of dissipative Landau-Zener 
tunneling very interesting. 

Also, the Landau-Zener tunneling is of interest in the 
context of quantum computing, because it describes the 
adiabatic flip of a qubit. During the process of the flip, 
the qubit is in a superposition of the pure "on" and "off" 
states and thus may be very sensitive to environmental 
noise. Thus, the problem of dissipative Landau-Zener 
tunneling may be important for the physical implemen- 
tation of the qubit. 

Previous works on the dissipative Landau-Zener tun- 
neling treated the problem either perturbatively or phe- 
nomenologically P,|7|-|l0|] . Perturbative treatment, how- 
ever, is insufficient for adiabatic limit, because it poses 
too strong limitations on coupling strength. Phenomeno- 
logical approach is satisfactory only for high tempera- 
tures. Also, previous works focused primarily on the case 
of Ohmic coupling which is critical for equilibrium tun- 
neling . We will argue that the dissipation is relevant 
for Landau-Zener problem when the coupling constant 
a has the same dimension as the velocity of levels rela- 
tive motion i/. (We call the latter bias in the following.) 
We call this coupling marginal. (A physical example of 
such coupling is an SNS junction with two Andreev states 
connected to resistor.) ft turns out that the tunnel- 
ing is blocked if the bias is small: v < a. For larger 
bias, we derive an instanton solution describing dissipa- 



tive Landau-Zener tunneling at low temperatures and 
find the tunneling probability: 



exp 



(1) 



where A is the energy gap, i' is the bias, and a is 
the coupling strength. In the high temperature limit 
(T ^ fi/S./2v) we derive master equation which takes 
into account both thermal noise and dissipation. This 
equation solves the problem for all temperatures larger 
than the crossover temperature, and can also incorporate 
non-equilibrium noise. 

Model. It is conventional to represent the two states 
of Landau-Zener model using the spin 1/2 basis. The 
Hamiltonian is 



H — via. 



He + H„ 



(2) 



where v is the bias, A is the energy gap, and (7x,z are 
Pauli matrices. The term He describes the coupling to 
the environment with Hamiltonian Henv 

If the coupling He is negligible, the model can be solved 
exactly 0|, and the non-adiabatic transition probability 
is w = exp(— ttA^/i/). In the adiabatic limit, ly <C A^, 
the "adiabatically frozen" eigenvalues are: 



e±(i)-±VPF+A^. 



(3) 



The characteristic time of tunneling is ro = A/i>. 

Following jlll l, we model the environment by a set of 
oscillators with coordinates Xi and momenta p^: 



Henv — ^ ^ 
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and couple the environment to the spin linearly: 

He=d,X{t)- X{t) =Y,l^X^{t) ■ 



(4) 



(5) 



We consider here only the diagonal coupling which is the 
main source of dephasing. The effect of environment de- 
pends only on the spectral function: 



1 



(6) 



Without the loss of generahty, in the hmit of smaU uj 
one can consider power-hke spectral functions: J{lo) = 
asUi''- To estimate the effect of dissipation, we look at 
the dimensionless ratio rj — Os/^i'Tq'^^), which gives the 
measure of the dissipation with respect to the bias v. In 
the adiabatic limit, when tq is large, the dissipation is 
weak (rj ^ 1) for s > — 1, and strong otherwise. The 
marginal situation arising when J — a/uj is a, subject of 
this work. 

Marginal coupling describes, e. g., single-channel SNS 
junction connected in series to the resistor. There 
are two Andreev states in the junction. The dynamics of 
these states is governed by the Hamiltonian 



SNS 



A / 0(^) / X 

cr^ + AoVl - T sm —^<7x , (7) 



where (j){t) is the superconducting phase difference across 
the junction, Aq is the superconducting gap in the leads, 
and r is the channel transmission coefficient. Let V is 
voltage drop across the whole circuit, and V]i{t) is the 
voltage drop across the resistor. Then, the phase differ- 
ence is: 



2eVt 



y I Vn{t)dt 



(8) 



In the vicinity of the point (f) = n, where the energy 
difference of the two states is minimal, one may expand 
HsNS in 0- If there were no resistor, the system then 
would be described by Landau-Zener model. To see what 
is the effect of the resistor, consider voltage fluctuations. 
According to quantum Nyquist formula, the voltage Vr 
fluctuates at zero temperature as (Vri^Vr^^ij) ^ Ruj, 
where R is the resistance, and lu is the frequency. Thus, 
phase fluctuations are: 



R 



(9) 



Comparing this correlation function to [Xi^X^J) , one 
may see that the circuit is described by Landau-Zener 
model with marginal dissipation. The parameters are 
identified as: 



sVAn 



A = AoVl~T 



ei?Ag 



(10) 



There are two regimes of tunneling for different val- 
ues of environment temperature T . For T ^ h/S./v the 
tunneling is thermally assisted, and the tunneling prob- 
ability obeys Arrhenius law: w ~ exp(— 2A/T). For low 
temperatures, the tunneling is due to non-adiabaticity, 
and the probability saturates at T ^ 0. We consider 
these regimes separately. 

Quantum regime. At low temperatures, we treat 
the nonequilibrium problem by the Keldysh technique 



||T2|,^. We introduce two time contours describing evo- 
lution of wave function and its complex conjugate. Stan- 
dard contours going along the real time axis are inappro- 
priate for the problem in question because the exponen- 
tially small adiabatic transition probability results from a 
delicate destructive interference of many oscillating con- 
tributions. We move Keldysh contours away from the 
real axis to avoid oscillating terms. Our choice of con- 
tours is similar to [Q, where it was used in the context 
of exciton autolocalization. To illustrate the use of con- 
tours, consider the tunneling for a = 0. The transition 
probability is w — exp(— ^o), with the action S given by 



5*0 = 2Im J (e+(t) - e^{t)) dt = 
ti 



(11) 



The integral is taken from any point ti on the real axis 
to itq, where the adiabatic energy e{t) has a square root 
singularity. The two adiabatic states correspond to the 
branches e±{t), i. e., to the two different sheets of the 
Riemann surface. 
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FIG. 1. (a) Keldysh contours, (b) Transformed contours. 
Dashed lines go on the other sheet of the Riemann surface. 
Wavy lines are branch cuts. 

The result (|ll]) may be represented using Keldysh con- 
tours shown on Fig||(a). In the Keldysh formalism, the 
forward contour always goes from t = —oo on the real 
axis, where the initial state is prepared, to t — +oo where 
the final state is measured. We draw the contour through 
the branching point t = itq. At this point the two en- 
ergies e+(t) and e_(i) coincide, hence the adiabaticity is 
violated, and the transition occurs. It means that we 
change the branch e_(t) a-nd continue the con- 

tour on another sheet of Riemann surface. Then, the 
contour goes back to the real axis, and then to < = +oo. 
The backward contour is a complex conjugate of the for- 
ward contour. In the adiabatic approximation, the seg- 
ments parallel to the real axis give no contribution to the 
imaginary part of the action, and thus one may connect 
vertical segments into a single contour C (see Fig|^(b)). 
In this representation, the action (ni]) is given by 



(12) 
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It turns out that the problem with a ^ has a solution 
of a similar structure: there is a square root singularity 
at some point t = ixa in the complex plane, and the 
branch is changed when the contour passes this point. 
In this case, the branching point iTa has to be found 
self-consistently. We derive the effective action on the 
contour C, look for an instanton solution and find 
from equations of motion. Let us introduce a field <I>(t) = 
lyt + X{t), so that the energy of the spin is a/A^ + ^^{t). 
Thus, the two-level system action is 



(13) 



(The Berry phase is neglected here since the "magnetic 
field" acting on the spin is always in the xz-plane). The 
dynamics of the environment is described by Senv + S\, 
where 



1 



2S' 



(14) 



and the term S\ enforces the constraint — X = u: 

Sx^-i£ dt\{t) (^{t) -V-Y. l^X^{t)^ 



(15) 



We integrate out $(t) and xiit) in the saddle-point ap- 
proximation and find the effective action in terms of the 
Lagrange multiplier X{t): 



Seff = iA <j> dt\J I - +w <j> dt\{t) 



— (J) 4) dt dr 

4tt 



{t~t') 



t\2 



(16) 



The solution describing Landau-Zener tunneling has 
different signs on different sheets of the Riemann surface: 
\{t — 0) = —\{t + 0). Taking a saddle point with respect 
to \{t), one arrives to the equation of motion: 



^ d X ia f \(t')dt' , , 

A— , . = = -t/+— / . (17) 



dt 



TT J {t-t'f 



For a = the solution is Ao(t) = \Jtq + 1'^. Inter- 
estingly, the integral in ( p^ ) with A = Ao(t) is con- 
stant for t € [— iTo;jro], and the r.h.s. of (^^ pre- 
serves its form. Thus, one may look for the solution 
of the form X{t) = y/r^TW. From Eq. (|l^ one finds 
Ta ~ A./{v — a). Substituting it into (|l6|), one obtains 
the action 



S 



7rA2 



V — a 



(18) 



and the tunneling probability exp(— S*) given by Eq. (Q). 



This solution gives an instanton for v > a. Otherwise, 
for V < a, there is no saddle point solution and the tun- 
neling is impossible. This happens because the friction 
force is larger than the bias, and the energy put into the 
system by the bias source, is fully dissipated into the en- 
vironment. Because of that, there is no level crossing, 
even at complex times. The tunneling time Tq, diverges 
ai V — a. Because of that, the dynamics far from level 
crossing point becomes essential. The linear approxima- 
tion on which Eq. ^ is based may break down. Also, for 
1/ ~ a, the tunneling time may be comparable to the time 
too when the system is prepared and measured. Thus, the 
Landau-Zener problem makes no sense for v too close to 
a. The behaviour of the system is determined by the de- 
tails of the dynamics far from the avoided crossing point. 

In the perturbative regime, a <^ the action ([l^) can 
be expanded in a: 



S — 5*0 ■ 



ttA^q; 



(19) 



This result contradicts to Ao et al., ||^, who found an 
ck-indepcndent tunneling probability. We believe that 
the result |^ is incorrect because of improper choice of 
shakeup force C,{t) that does not describe tunneling. 

The above solution shows that at T = the envi- 
ronment effect is purely dissipative, and the coupling to 
the environment reduces the transition probability. This 
happens because zero-point fluctuations of the environ- 
ment do not lead to a real transition. At finite tem- 
peratures, the environment may transfer the energy to 
the two-level system, thus increasing the transition prob- 
ability. If the temperature is low, T <^ Tq — Ti/Ta^ 
the transition still has mostly quantum nature. In this 
limit, the finite temperature can be taken into account 
by requiring the instanton to be periodic in imaginary 
time, \{t) = \{t + iP), with the period P = 1/T. 
Thus, one has to consider a periodic system of cuts, 
t e [—IT -|- iPn, IT -\- i(3n\. Note that the nonlocal term in 
( |T6| ) gives rise to interaction between instantons on dif- 
ferent cuts. If /3 ^ Tq,, the interaction is weak and can 
be treated perturbatively. The correction to the action 
per period is: 



5S 



1 



n#0 ^ 



After computing the sum over n, one obtains: 



S{T) = S'(T = 0) 1 



6 {v — a) J 



(20) 



(21) 



and the tunneling probability is w ~ exp(— 5(r)), as 
before. Eq. ( pl[ ) is applicable when the correction is 
small and thus does not lead to non-monotonous T- 
dependence. 

Classical regime. If the temperature of the environ- 
ment is high (T ^ h/Ta), the tunneling is thermally 
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assisted, and the effect of environment can be divided 
into slow regular motion with frequencies lo ^ 1/tq and 
fast Langevin noise with lo ^ 1/tq. (Only noise with 
Lj > 2A contributes to the transition probability, and 
this separation is valid in the adiabatic limit.) Then the 
Hamiltonian is: 



H = F{t) + Ao-^ + u{t)a^ 



(22) 



where F{i) = vt + {X{t)) is a regular part, and u{t) — 
X{t) — {X{t)) is a fluctuating part. Since the first two 
terms are slow functions of time, one may consider them 
in the adiabatic approximation, and treat the noise u{t) 
as a perturbation causing non-adiabatic transition. This 
implies that the noise contribution to the transition am- 
plitude is larger than the non-adiabatic correction. The 
eigenstates of the frozen Hamiltonian are: 



cos ■ 



V'-(i) 



(23) 



where ta.n 6{t) — A/F{t) . In this basis the Hamiltonian 
is: 

H = -e{t)d;^ +u{t){cos0{t)a^ +sm9{t)a^) , (24) 

where e{t) ~ A/ smd{t) is the adiabatic energy. The per- 
turbation theory with respect to u{t) gives the transition 
rates 



^ib,m — >n — 



smd{t) {n\u{t)\m) dt 



(25) 



-oo — oo 



(26) 

This integral contains a fast oscillating function, and is 
determined by the region where \t — t'\ ~ l/e(t) ^ tq. 
Therefore, one may approximate in the prefactor t t' , 
integrate over t~t' and find the Landau-Zener transition 
rates: 



oo 

u;± = y K{uj = ±2e{t)) sin^ 9{t) dt 



(27) 



is the noise spectrum. In equilibrium, it can be expressed 
in terms of the spectral function J(a;) and the Bose dis- 
tribution function N{lu): 



K{lo)^J{\lo\){N{\lo\)+9{u;)) . 



(29) 



The result (^7|) is also true for non-equilibrium noise with 
the proper choice of K{lu). 

To find the still unknown function 0{t), consider the 
force exerted by the rotating spin on the environment, 
f{t) — {<7z{t)) — ±cos9{t). The response to that force 
is, in Fourier representation, (X^) = —iJ{bj)f^. Then, 
one may write an equation for 9{t): 

F(t) = iy+{X{t)) ^vTacoseit) . (30) 

Since F{t) — A/tan0(t), one has: 

Ae = sm^ 6 {t){iyT a cose (t)) . (31) 

For small bias, v < a, there is no Landau-Zener tunnel- 
ing, and situation is similar to quantum limit. If the bias 
v > a, the transition probability is: 



w± - / r±i0)de ; r±{0) 



V cos 9 



(32) 



between the states |m) and \n) of the environment. In 
Eq. (HI), \{t) — e{t), and the subscript (±) denotes 
initial spin state. Tracing out the environment, one finds 
the Landau-Zener transition rates 

oo oo 

w±^ [ dt [ di'e=^'(^(*)-^(*'» sinei(i) sin6l(t')(^t(0"(^')) • 



Note that the small bias v appears in the denominator. 
Because of that, the probability computed from ( |3^ ) can 
be quite large, and the perturbation theory on which this 
equation is based, appears to break down. This happens 
because the motion at small bias is slow, and a small 
noise acts on the system for a long time, giving rise to a 
large tunneling probability. 

To solve the problem, note that according to (p7|), there 
is no interference between transition amplitudes at differ- 
ent times t. This happens because the transitions are due 
to random noise, and the amplitudes of separate transi- 
tions have random phase. In this situation, one may use 
master equation, with the transition rate per unit time 
from (p7|). Since 9(t) depends on the history of the sys- 
tem, it is convenient to use 9 as an independent variable 
instead of t. Master equation, written in the variable 9, 
has the form: 



dP+{9) _ dP^{9) 



d9 



d9 



^r+{9)P+{9)+r_P_{9) , (33) 



Here 



where P± (9) are occupancies of two adiabatic states. The 
initial condition is F+(0) = 0, P_(0) = 1. 

Eq. (^^ can be solved in a general form. However, 
for simplicity, we consider the two cases: (i) T <C A and 
(ii) T 3> A. For T ^ A, the transition rates r_|_ and r_ 
are equal. The noise correlator is K{uj) = aT juP' . The 
calculation gives for w = P+ (6* = tt) : 



(u{t)u{0)) dt 



(28) 



1 



(1 



-2R\ 



R = 



^{9)d9 . 



(34) 
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Computing the integral, one obtains: 



„ ttT I V 



(35) 
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In the hmit T ^ oo, Eq. (|^ predicts w —f 1/2, the 
result of a phenomenological approach. Note that (^) 
is true only when R <^ A^/:^, otherwise non-adiabaticity 
is again important, and the exponential correction has 
different structure |^,||. 

We consider now the case of intermediate tempera- 
tures: fi/Ta < T < A. The transition l^;-) 
requires energy transfer 2A, while the inverse transition 
does not, and the transition rates r_(_ and r_ are not 
equal. Therefore, the transition from the lower state to 
the upper one occurs at 9 = 7r/2, when the energy differ- 
ence is minimal. This transition can be considered as an 
instant excitation occurring at ~ 7r/2 during the time 
tea ~ V AT/i/. It is followed by a decay from the upper 
state. Thus, the transition probability is the product of 
the probability of excitation at 6 = Tr/2 and the prob- 
ability that the excited state survives for 7r/2 < 9 < tt. 
From the master equation ( |3^ ) one finds the excitation 
probability: 



p — 



r_{9)d9^-xl^e'^^/^ 



(36) 



and the survival probability: 
/ . 

P, = exp - y r+{9)d9 
Finally, the transition probability is: 



(37) 




-2A/T 



(38) 



Note that the dependence of this expression on a is non- 
monotonous. This is due to the fact that the coupling to 
the environment provides both external noise and dissi- 
pation. The transition probability reaches the maximal 
value at a = 1v j'i. 

Conclusion. We studied dissipative Landau-Zener tun- 
neling marginally coupled to environment in both quan- 
tum and classical limits. At T = the effect of environ- 
ment is purely dissipative, whereas at high temperatures 
thermal fluctuations in the environment increase the 
transition probability, leading to non-monotonous de- 
pendence of the transition rates on the coupling strength 
a. At v < a the tunneling is blocked. Thus, the inter- 
play between fluctuations and dissipation is manifest in 
the marginal coupling model. 
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